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ABSTRACT 

Dynamical studies of local elliptical galaxies and the Fundamental Plane point to a 
strong dependence of the total mass-to-light ratio on luminosity with a relation of 
the form M/L oc L 1 . The "tilt" 7 may be caused by various factors, including stellar 
population properties (metallicity, age and star formation history), IMF, rotational 
support, luminosity profile non-homology and dark matter (DM) fraction. We eval- 
uate the impa ct of all these facto r s usin g a large uniform dataset of local early-type 
galaxies from Prugniel & Simicn (1996). We take particular care in estimating the 
stellar masses, using a general star formation history, and comparing different popu- 
lation synthesis models. We find that the stellar M/L contributes little to the tilt. We 
estimate the total M/L using simple Jeans dynamical models, and find that adopting 
accurate luminosity profiles is important but does not remove the need for an addi- 
tional tilt component, which we ascribe to DM. We survey trends of the DM fraction 
within one effective radius, finding it to be roughly constant for galaxies fainter than 
Mb ~ —20.5, and increasing with luminosity for the brighter galaxies; we detect no 
significant differences among SOs and fast- and slow-rotating ellipticals. We construct 
simplified cosmological mass models and find general consistency, where the DM tran- 
sition point is caused by a change in the relation between luminosity and effective 
radius. A more refined model with varying galaxy star formation efficiency suggests a 
transition from total mass profiles (including DM) of faint galaxies distributed simi- 
larly to the light, to near-isothermal profiles for the bright galaxies. These conclusions 
are sensitive to various systematic uncertainties which we investigate in detail, but are 
consistent with the results of dynamics studies at larger radii. 

Key words: dark matter - galaxies : evolution - galaxies : galaxies : general - 
galaxies : elliptical and lenticular, cD. 



1 INTRODUCTION 

Early-type galaxies (ETGs) are the most massive stellar sys- 
tems in the Universe, containing much of the cosmic budget 
of visible and dark matter (DM). They include elliptical (E) 
and lenticular (SO) galaxies and form a nearly uniform class 
of objects: usually red, old and with only traces of cold gas 
and active star formation. The striking regularities in their 
properties include strong correlations between size (e.g. the 
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effective radius, R e g) and the surface brightness therein (I e s; 
iKormendvl Il977h . and between kinematic s (the central ve- 
locity dispersion 00) and luminosity (L; IFaber fc Jacksonl 
11979 . hereafter FJ). 

The two relati ons above merge into the so-called Funda- 
ment al Plane (TP: iDiorgovski fc Davisl ll987: Dr essier et all 
Il987i ). i.e. a relation between the (logarithm of) o~o, R c a 
and /off of ETGs. The FP can be interpreted in terms of 
the virial theorem of relaxed systems, according to which 
2T + U — where U is the potential energy and T the ki- 
netic energy. This can be re-written in terms of observed 
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quantities as approximately L oc a 2 R c g. However, the FP 
is found observationally to be L oc a v R c ft a with a 7^ 1 
and 77 7^ 2, i.e. a different orientation of the plane in the 
space of the logarithmic quantities with respect to the virial 
prediction. This tilt of the FP provides insight for the for- 
mation and structure of ETGs, and can be interpreted as 
a variation of the to tal mass-to-light ratio (M/L) with L 
jDressler et al.ll 19871) with the simplest parameterization as 
a power law, M/L oc L 1 . The slope, 7, of this relation 
could be driven by one or more different factors: a vari- 
ation in stellar M/L (due to metallicity or age gradient 
or change in IMF), a variation in the D M content, non- 
homology, rotational support, etc. (see e.g . iBusarello et al.l 
Il997l ; ID'Onofrio et~aT] I2OO6I ; iGravesI [iooi l . It is of consid- 
erable importance to disentangle these factors using high- 
quality data at low redshift, in order to use the FP as a 
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Many studie s over the years have attempted to decode 
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see however ljun fc Imll2008l ). with the major contributor 
yet to be firmly identified-which would have r amifica- 
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An additional complication is that the tilted FP may 
not be flat, with claims made for curvature llZaritskv et al.l 
120061 ; ID'Onofrio et~aH |200S| ; iHvde fc Bernardill2008bl ). and 
projections of the FP showing a bend at a characteris- 
tic magnitude of Mb ~ —20.5 (see Section [2}. This tran- 
sition concords with more general findings of a discon- 
tinuity in ETG properties at a similar luminosity (e.g. 
iNapolitano et al.ll2008l ; ICoccato et al.ll2009l ). 

Of particular interest is the FP contribution from 
the central DM content, since connecting luminous galax- 
ies to their DM haloes is one of the key ingredients in 
modern recipes for galaxy assembly. The DM component 
of the FP tilt is currently a bone of co ntention, with 
findin gs alternatively of negligible impact (|Gerhard et al.l 
l200ll. hereafter G+01 ; T+0 4 ) and of primary im portance 
ijPadmanabhan et al.l 120041 ; iLintott et all l2006l ; C+06). 
More generally, several lines of evidence point to a small 
but sig nificant fraction of DM (/dm) inside 1 R e g (e.g. 
G+01; iGavazzi et ail l2007l : IWeiimans et al.1 I2008I ). but 
the galaxies probed are typically very bright Es with 



(To > 200 km s . As data have become available on 
more "ordinary" ETGs, there are suggestions that their 
DM properties vary strongly with luminosity, with per- 
haps even a dichotomy following the classic division of 
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The purpose of this paper is to reexamine all the plausi- 
ble factors that could contribute to the FP tilt. We will not 
be directly studying the FP, but rather will survey system- 
atic trends in the central properties of ETGs which factor 
into the apparent M/L variation, 7. Our study springs from 
the classic data-set of central photometry and kinematics of 
PS96, which is one of the largest such catalogues of local 
galaxies covering a wide range of luminosities. Our analysis 
consists of two parts, wherein we independently determine 
the stellar component of the M/L by stellar populations 
models, and the total M/L by dynamical models. We take 
particular care in considering realistic star formation histo- 
ries, stellar populations modelling systematics, and dynam- 
ical contributions from DM. 

We briefly describe the galaxy sample in Section [2] and 
analyze the stellar populations in Section [3] We determine 
the dynamical masses in Section [4] and infer the DM frac- 
tions in Section [5] Section [6] considers some implications for 
galaxy formation, and Section [7] draws conclusions and con- 
siders future prospects. 



2 SAMPLE 

Our data-set of local ETGs is drawn from PS9(£]. This is 
currently one of the largest homogeneous samples of local 
ETGs available in the literature containing both photom- 
etry and kinematics of the galaxy central regions, and the 
only one including information on the peak rotation velocity 
(Knax). As discussed in PS96, the colours (extinction- and 
K-corrected) are measured within 1 _R e ff , the central velocity 
dispersions 00 are recovered from long slit spectr43 and V max 
is defined as the quadratic sum of the maximum rotation on 
the major and minor axes. Since we are interested in fitting 
spectral energy distributions (SEDs), we select galaxies with 
at least two measured colours (most of the selected galaxies 
have four colours jfl Selecting also for galaxies brighter than 
Mb = —16, we recover ~ 400 galaxies among which, follow- 
ing the PS96 classification, ~ 55% are bona fide Es (their 
subsample 1), « 30% are type SO and Sa (their subsample 



1 Downloadable at http://vizier.cfa.harvard.edu/viz 
-bin/VizieR?-source=J/A+A/309/749. 

2 The authors do not report detailed information about their 
measurement set-up but, as extensively adopted in literature 
analyses of this dataset, we will interpret <tq as the luminosity- 
weighted velocity dispersion within a circular aperture of radius 
Be«/8. 

3 Apparent total magni tudes are on average slightly brighter than 
those in RC3 catalog (de Vaucoulcu rs et al1ll99ll) by —0.05 ± 
0.1 mag, while the differences of B — V and U — B colours with 
those in RC3 are 0.00 ± 0.03 and -0.02 ± 0.03. 
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Figure 1. Correlations among photometric and kinematical properties of Es in our galaxy sample. The small points show individual 
galaxies, while the large points with error bars show binned averages. The solid lines show linear fits, carried out separately for the bright 
and faint subsamples (divided at Mb ~ —20.5, i.e. log Lb ~ 10.4 Lq). Left Panel. FJ relation. Right Panel. Size-luminosity relation. 



5), and the remaining ~ 15% are dusty objects, interacting 
galaxies, dwarf spheroidals, compact, dwarf, low-luminosity 
and peculiar ellipticals, etc. (subsamples 2, 3, 4 and 6). For 
the main purposes of this paper, we will use subsample 1 
(hereafter "Es") and 5 (hereafter "SOs"), thus incorporating 
335 galaxies, or « 85% of the PS96 sample. 

In all the following, we use a cosmological 
model with ({l m ,£l\,h) = (0.3,0.7,0.7), where 
h = ffo/lOOkms- 1 Mpc" 1 (|Spergel et all I gOOSl ). corre- 
sponding to a universe age of t un iv = 13.5 Gyo 

E galaxies populate a restricted region in the colour- 
magnitude diagram, the so-called "red sequence", with a 
colour range of B — V ~ 0.9 — 1. SOs span a wider range of 
colours (i.e., B — V ~ 0.7 — 1) and are fainter than Es on 
average. The two subsamples follow similar FJ relations L oc 



ot/ 1 , including a characteristic magnitude {Mb 



-20.5) 



where the relation clearly changes its slope (see left panel of 
Fig.Qjj. For the E sample, tj = 2.9 ±0.5 and 5.6 ±1.2 for the 
faint and bright galaxies, respectively. 

This FJ "dichotomy" has been reporte d elsewhere 
l|Matkovic fc Guzman! [20051 ; iForbes et aHl2008l ). and seems 
related to systematic changes seen i n other optical 
prop e rties, e.g. the Kormendy relati on dCapaccioli et all 



l2007l ; lHvde fc Bernardill2008al ). The latter relation is illus- 
trated for our sample Es in the right panel of Fig. [T] where 
the faint and bright galaxies have fitted slopes of 0.85 ±0.07 
and 0.43 ± 0.11, respectively. 

We will see that a characteristic luminosity scale is also 
found to characterize other correlations of ETG parameters. 



3 STELLAR MASS-TO-LIGHT RATIO 

One of the key aspects of our analysis is determining each 
sample galaxy's stellar M/L, T*, which we do by fitting 
model SEDs to the observed galaxy colours. Although pho- 
tometric modelling may seem less powerful than detailed 
spectroscopic fits, most spectroscopic samples are restricted 
to the very central regions of galaxies and may be very bi- 
assed indicators of the stellar populations on scales of ~ J? e ft . 

In Section [371] we describe the modelling procedure and 
present the recovered stellar populations properties for the 
sample galaxies. We report the implications for the size-mass 
relation in Section f3. 21 and for trends in T* in Section f3. 31 
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4 The distance scale is critical for normalizing the luminosities 
and M/Ls. The distance moduli (m — M) from PS96, rescaled 
to h = 0.7, are on average lower those reported in the RC3 cat- 
alog llde Vaucouleurs et al.|[l99lh by — O.II^q'jq mag (uncertain- 
ties are 25 th and 75 th percentile s), while agreeing closely with es- 
timates from lTonrv et all ll200l|) (shifted by —0.0 6 mag to correct 
to the Cepheid distance scale; I Jensen et al.l2003h which differs by 
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+0.18 



mag. 



We create a set of synthetic st e llar s pectra using the 
prescription of iBruzual fc Chariot! (|2003l . hereafter BC03), 
which encompasses a wid e range o f initial metallicit ies and 
ages. A lSalpeterl (|l955l ) or lChabrierl (|200ll . l2002l . l2003l ) initial 
mass function (IMF) is assumed, with initial masses m in the 
range 0.1 — 100. The two IMFs do not influence the colours, 
but basically affect the T* estimates, which are ~ 1.8 times 
larger with a Salpeter IMF than with a Chabrier IMF. 

To generate a more general and realistic star forma- 
tion history (SFH), we convolve the BC03 "single burst" 
models with an exponentially-decaying star formation (SF) 
rate with time t, oc e~ l ^ T , where r is a characteristic time 
scale. The choice of BC03 is dictated mostly by its versatility 
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Figure 2. Monte Carlo simulation of the stellar populations fit- 
ting procedure, where the estimated stellar M/L values T es t are 
plotted against the intrinsic model values T s ; m . The solid line 
is the one-to-one relation. The input value of T» is reproduced 
with no systematic bias, and the scatter of the recovered values 
is ~ 10%. 



and ability to span the stellar parameter space (metallicities 
and ages) but it is not the only prescription available on 
the market. We test for the presence of any modelling sys- 
temat ics by mainly checkin g two different popula r prescrip- 
tions, iBell fc de Jonj (|200ll , hereafter BdJOl) and lMarastonl 
|2005l . hereafter M05) in Appendix [Aj 

For each galaxy, we fit synthetic spectra to the observed 
colours (U — B, B — V, V—R and V — I, after convolving the 
spectra with the appropriate filter bandpass functions), al- 
lowing us to estimate the age (t ga i), metallicity (Z), r, T» t s 
and hence the stellar mass, M*= T,^ x Lb (hereafter we 
will always quote luminosity and M/L values in the B-band, 
even if not specified). In detail, we build a set of synthetic 
colours with Z € (0.008, 0.02, 0.050, r <E (0.1 - 5) Gyr and 
igai up to tuni-tHl The fitting procedure consists of generating 
100 Monte Carlo realizations of the observed galaxy colour 
sets assuming Gaussian errors of 0.05 mag per colour, and 
minimizing a x 2 statistic between the modelled and observed 
colours for each realization. The overall best-fit model pa- 
rameters and their uncertainties are defined as the median 
and scatter of these 100 best fit 

Our synthetic modelling procedure is more general 
than the extensively used "simple stellar population" (SSP) 
model where a galaxy is approximated as experiencing a 
single burst of star formation (i.e. r = 0; iTrager et al.ll2000l: 
Maraston 2005; see also Appendix|X]|- Instead we leave t ga i, 



5 Lower metallicities would have produced t ga j larger than the 
age of the universe in our assumed cosmology (~ 13.5 Gyr). 

6 We allowed ~ 10% scatter around this value in order to account 
for some intrinsic uncertainty in the age estimates, and thus some 
of the estimated ages might be slightly larger than 13.5 Gyr. 

7 Typical la uncertainties on the estimated T„ are ~ 10 — 20%. 
In general in this paper, we will fit medians rather than means in 
order to be more robust to outliers. 



r and Z all as free parameters in order to better repre- 
sent the wide variety of SFH expected both observationally 
and t heoretically (see e.g. lDe Lucia et al . 2006; Noesk e~et al.l 
2007). The allowed ranges in the parameters will then be 
larger than in the more simplified SSP case because of the 
well known degeneracies among them (|Gavazzi et al ] |2002l ; 
BC03). 

To test the reliability of our modelling technique and 
the intrinsic parameter scatter, and to check for the pres- 
ence of spuriously-generated correlations, we run a suite of 
Monte Carlo simulations. We extract 100 simulated galaxy 
spectra from our BC03 SED libraries with random t ga i, Z 
and r (i.e. with no correlation among these parameters), and 
apply our generalized fit procedure — comparing the recov- 
ered parameters with the input model values. We find that 
T„ is recovered well, with a scatter of ~ 10% (see Fig. [2}. 
Similar consistency is found for t ga i, Z and r, which have 
on average larger scatter: ~ 20%, ~ 30% and ~ 30% respec- 
tively. We check for spuriou s correlations using a Spearman 
rank test (|Press et al.lll992T ). finding that r vs t ga i show no 
correlation at the 95% confidence level, but that t ga i and Z 
are weakly correlated, which is a common effect in stellar 
populations analyses, However, this t sa ,\-Z degeneracy does 
not affect the T* inference, which is our primary concern. 

After fitting the real data for the complete sample of 
335 galaxies, we show some relations between model param- 
eters and other observed galaxy quantities in Fig. [3] From 
this Figure it is evident that the metallicity is generally so- 
lar or super-solar [Z ^ 0.02) and on average only weakly 
dependent on other properties such as luminosity - which 
is fortunate since the BC03 stellar libraries include only a 
few reference values for Z. More striking are the strong cor- 
relations involving t ga i and t, such that the brighter, more 
massive galaxies formed their stars on average on shorter 
timescales than the fainter, less massive galaxies, while the 
younger galaxies also had shorter SF timescales (this does 
not mean that the brighter galaxies are younger, and if the 
SOs are included, the opposite is clearly true). 

Similar findings on the SF Hs of ETGs have 
been found i n ob servational (e.g. iGavazzi et alj 20021; 
Thomas et al. 20051) and theoretical (e.g. |Pe Lucia et al.l 



20061 ; iRomeo et al.l [2008) analyses. We will consider this 



subjec t in detail in a subsequent paper (jNapolitano et all 
2009b, in prep.), and for now summarize some basic pa- 
rameters. The bright ETGs (Mb ^ —20.5) have a median 
t ~ 0.5 Gyr, while the faint ETGs have r ~ 1 Gyr. The Es 
have similar properties, while the SOs have on average more 
protracted SFHs (r ~ 1—1.5 Gyr). The median T* for the 
ETGs is 6.9 ± 2 T Q (3.8 ±1.1 Tq) for a Salpeter (Chabrier) 
IMF, where the quoted errors are the 1 a scatter. The SOs 
have only slightly smaller median T» than the Es: 6.2 Tq 
(3.4 Tq) and 7.1 Tq (3.9 Tq), respectively. However, the two 
samples differ more strongly in the distribution of T» , where 
the Es have a fairly symmetric distribution about the mean, 
while the SOs have a pronounced tail to low T* such that 
the mean value of 5.8 Tq (3.2 Tq) differs from the median 
estimate. 

3.2 Size-mass relations 

An indirect way to test our derived T, values is to check 
how the implied scaling relation between size and stellar 
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Figure 3. Stellar populations properties of E sample galaxies: mctallicity Z, SF timescale t, stellar age i ga i and luminosity Lb- In each 
panel the data are binned in intervals, with the median and ±25% scatter shown. The fraction of galaxies in each bin is reported in the 
first two panels. For SOs, the trends with metallicity are identical to the Es, while the ones with r are weaker. 



mass compares to previously established results. Expressing 
the size- luminosity relation as R c g oc L aL , we find a slope 
for the Es of = 0.70 ±0.06, which is s lightly steeper than 
some literature findings of 0.54 — 0.63 dPahre et alj [l998a , 
iBernardi et alj|2003l . [Mamon fc Lokasll2005r i. However, this 
slope is very sensitive to the range of luminosities fitted, 
since we see a difference between the "faint" and "bright" 
subsamples (see Section [2] and right panel of Fig. [T}. 

For the size-mass relation R c s oc M+ M , we ob- 
tain (x m — 0.65 ± 0.05 overall (Fig. Q, consistent with 
pre vious estimates for low-redshift galaxies of typically ~ 
0.6 jBernardi et alJliool ; IShen et al.|[2003l ; iMamon fc Lokasl 
12003 ? N+05). This correlation also bends at a characteris- 
tic mass scale of M» ~ 1Q 1L1 M S , with a M = 0.36 ± 0.13 
and 0.73±0.12 for th e faint and bright galaxies, respectively 
(cf. IShen et al.ll2003T ). We find the SOs to have on average 
smaller R e ff than Es of the same mass, and similar aM at 
high masses, but flattening to ~ at low masses. 



3.3 Luminosity dependence of T* 

The central regions of ETGs are probably dominated by 
the stellar mass, so it is of critical importance to ascertain 
the fraction of the FP tilt that is connected to the stellar 
population properties. We focus on the relation T* oc Lb 1 * 
in log-log space, fitting to weighted medians of binned data 
values, with results that are stable to changes in the binning 
(Fig. [5}. For the overall ETG sample, we find a slope of 
7* = 0.06 ±0.01; the Es have 7, = 0.02 ±0.01, and SOs have 
7* = 0.17±0.03, although this steeper slope is driven by the 
very faintest galaxies (Mb > —19). 

Before taking these results at face value, we investigate 
possible dependencies on modelling systematics. Changing 
the IMF from Salpeter to Chabrier does not affect 7*, 
but only the overall normalization of T». Adopting a sim- 
plified model with i ga i as the only free parameter, with 
r ~ 0.75 Gyr and Z ~ Z@ fixed to the median values for 
the whole sample (see Section f3 . 1 f) . the T* steepens for the 
faint galaxies and flattens for the bright ones, with an over- 



6 Tortora et al. 



1.25 



0.75 




-18 -19 -20 -21 -22 -23 



0.25 



-0.25 



Figure 4. Correlation between effective radius and stellar mass 
in E galaxies. The symbols are as in Fig. [T] The dashed line is 
the best-fit relation over the full range of the data, while the solid 
line is the best-fit to the two mass regimes. 



all result of 7, ~ 0.16 ( Appendix |A| . An even more simpli- 
fied SSP model with r = yields about the same 7* ~ 0.18. 
Thus, we see that allowing for the variations of SF timescales 
(and metallicity) within the sample is critical to accurately 
deriving the T* trends with luminosity. An additional com- 
plication, only partially addressed by our model's protracted 
exponential SFH, is multiple bursts of SF in a single galaxy 
which, if corrected for, would probably flatten the slope even 
further (cf. section 4.7 of C+06). 

We next consider alternative stellar populations basis 
models. As detailed in Appendix [X] adopting the BJ01 or 
M05 models with the same assumptions would imply shal- 
lower and steeper T* slopes, respectively. This turns out to 
be the dominant source of uncertainty in the analysis, al- 
though the uncertainty is largest for the faintest galaxies 
(Mb J> — 20.5), and the T, estimates more stable for the 
brighter galaxie^f]. 

How do our results compare to previous work? As men- 
tioned in Section [1] most recent studies agree that 7* is a 
relatively small contributor to the total 7 of the FP. E.g. 
PS96 inferred from the same data set that 7* ~ 0.1, and 
T+04 found using very different data and techniques that 
7* = 0.07 ± 0.01 in the B-banc(f|. This general consistency 
is encouraging, but it should be kept in mind that as just 
discussed, there remain significant uncertainties in the mod- 
elling. 



8 Including all the systematics, the final slope becomes 7* = 
0.06 ± 0.01+° ll. 

9 For the modelling, PS96 used fitting relations linking single 
colours and line-strength index Mg2 to the magnitude. T+04 ap- 
plied the PEGASE prescription l|Fioc fc Rocca-VoImerage![l997l) 
to the Sloan Digital Sky Survey (SDSS) Early Data Release. T+04 
also found 7* = 0.02 ± 0.01 in the /f-band using 2MASS data 
llBell et al.ll2003ft , but one would expect 7* to depend on band- 
pass because of the changing contributions from the mix of stellar 
populations. 
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Figure 5. Stellar M/L as a function of B-band luminosity, for 
Es (top panel) and SOs (bottom panel). Grey points represent in- 
dividual galaxies, while points with error bars are binned median 
values and scatter. Solid lines are the best fit linear relations in 
log-log space. A Salpeter IMF is assumed. 



The foregoing conclusions are based upon a universal 
IMF, but there are suspici ons that the IMF may vary wit h 
time or environment (e .g. Ivan Dokkumll200i ; |Pavell2008f i. 
iRenzini fc Ciottil (|l993h pointed out that a variation in the 
IMF with luminosity could easily account for the FP tilt. 
Here we illustrate this point again with a simple toy model, 
wherein the faint galaxies have a Chabrier IMF, and the 
bright ones a Salpeter IMF (a more realistic scenario would 
have the IMF changing smoothly with luminosity). The im- 
plied T» slope would be 7* ~ 0.3 (Fig. [B), which as we will 
see would be enough to explain the FP tilt with no further 
ingredients (e.g. no DM). 

We lastly examine the correlation of T* with velocity 
dispersion: T» oc ao~"° . The fitted slope is "/ a ~ 0.2 — 0.4 for 
both Es and SOs (Fig. [TJl . The steepness of this trend relative 
to 7» suggests that the stellar populations of galaxies are 
more strongly linked to their dynamical masses than to their 
luminosities. This issue will be considered in more detail in 
the following Sections, an d as part of a separate analysis in 
iNapolitano et al.l (|2009bl ). 



4 DYNAMICAL MASS 

Besides T„, the other fundamental quantity we want to de- 
termine is the total dynamical M/L, Td yn - The usual way 
the dynamical mass is calculated in FP studies is with the 
virial relation 



M 



Kalr 



(1) 



G ' 

where G is the gravitational constant and if is a 
pressure correction coefficient (or virial coefficient; e.g. 
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Figure 6. Stellar M/L for ETGs as a function of _B-band lumi- 
nosity, assuming an arbitrary change in the IMF at intermediate 
luminosity log(L/LQ) = 10.4. The faint galaxies (grey points) 
have a Chabrier IMF and the bright galaxies (black points) a 
Salpeter IMF. The blue solid line is the best fit, while the dashed 
one is the relation with a constant Salpeter IMF. 
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Figure 7. Stellar M/L ratio as a function of velocity dispersion 
(TO, using Salpeter IMF for Es (top panel) and SOs (bottom panel). 
See Fig. [5] for colour code and other details. 



iPadmanabhan et al.|[20o3 ; lEke et"al]|2004 C+06). The K 
factor subsumes a host of issues such as the aperture wherein 
(Jo is measured, and variations or "non-homologies" in the 
viewing angles, orbital structures, luminosity profile, DM 
distribution, etc. 

Rather than adopting some approximate value or func- 
tion for K , we will here directly model the non-homologies as 
much as possible for every individual galaxy in the sample, 
using the Jeans equations to estimate Td yn within 1 R e s. 



We outline the modelling methods in Section 14.11 present 
results based on luminosity profile homology in Section \4. 21 
and on more general profiles in Section [4.31 



4.1 Dynamical methods 

The basic approach of our dynamical models is to take the 
observed luminosity profile for each galaxy, along with a pa- 
rameterized mass model, and solve for the projected velocity 
dispersion ao within a central aperture. The mass model pa- 
rameters are then optimized to match the observed value for 
cr . 

In detail, the steps are the following: 

(i) We paramete rize the luminos i ty pr ofile j*{r) by ei- 
ther a ( deproj e cted)lde Vauco ulcurs ( 19 4^) profile or a more 
gener al ISersid |l968) model ( Caon. Capaccioli fc D'Onofriol 
Il993h -which fully takes into account any non-homologies 
in the stellar density distributions. The functional form for 
j*(r) is specified in Appendix B of PS97. 

(ii) We adopt a simplified form for the total cumulative 
dynamical mass profile M(r) which is either a constant- 
M/L profile M(r) = To L(r) (including th e cases where 
DM is missing or has a c ored distribution; see lBurkerj[l995l ; 
iNapolitano et akl |200 9a) , or a singular isothermal sphere 
(SIS), where M(r) oc 0"sis r - The latter choice is motivated 
by evidence from strong gravitational lens ing for near-SIS 
profiles in the central r egions of ETGs (e.g. lKocha nek 1991; 
iKoopmans et al1l2006l ). These two alternatives bracket the 
plausible range of mass profiles. 

(iii) We solve the Jeans equation: 

d(j*a 2 r ) , - fl(r) . 2 . GM(r) 



dr 



(2) 



where j3 = 1 — a 2 /of is the anisotropy. Th is model assumes 
spher ical symmetry and no rotation (cf. iMamon fc Lokasl 
I2005T ). For simplicity we also assume isotropy ((3 = 0), in 
which case the Jeans Eq. @ can be transformed to: 



a r (r) = 



. GM . 



(3) 



(iv) We project Eq. (J3J) to obtain the line-of-sight velocity 
dispersion: 



where 
I(R) = 2 



2 f°° j*<Tr rdr 



3* r 



-R 2 



zdr 



(4) 



(5) 



is the projected density profile. 

(v) We integrate <ti os within a fixed aperture R e s/8 to 
obtain the aperture velocity dispersion, <ta p using the Equa- 
tion: 



L(R) J 



Ree/8 



2nSI(S) a? os (S) dS 



(6) 



where L(R) — J Q 2nSI(S) dS is the luminosity within the 
projected radius _flF°l. 



We check that using an aperture of R e ff/10 would leave the 
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(vi) We fit the model cta p to the observed oo and iter- 
ate the preceding steps, varying the free parameters in Eq. 
([3]) (i.e. ersis or To). The resulting best-fit mass profile then 
provides the total spherical mass-to-light ratio within an ef- 
fective radius Tdyn(^eff) (which is coincident with To in the 
case of the constant-M/L model). 

This procedure does not take into account certain fac- 
tors that could in principle alter the final mass estimates. 
Firstly, the mass model does not include a central black hole, 
but we calculate the effect to be neglig 

ibkB 

More impor- 
tantly, real galaxies are neither spherical nor isotropic in gen- 
eral. We will check the impact of these simplifications later, 
but here begin with a first-order correction to the isotropic 
results. 

Detailed dynamical models of nearby galaxies have 
shown that their central stellar parts are close to isotropic af- 
ter subtracting the rot ational component (e.g. G+01; C+06; 
ICappellari et al.ll2007h . The observed cta p does incorporate 
both the projected rotation and dispersion components of 
the specific kinetic energy (<r Ap = « r 2 ms = v 2 + a 2 ), and the 
Jeans equations could in principle be reformulated along 
these lines. However, for many galaxies the rotation is so 
dominant that it is preferable to include it as an additional, 
separate term, which would require additional assumptions 
about the rotation field of each galaxy, and would best entail 
a non-spherical treatment anyway - all of which is beyond 
the scope of the current paper. 

Here we adopt a heuristic correction to the observed 
dispersion in order to approximately account for rotational 
effects. Following PS94, we parameterize the corrected er Ap 
by a' Ap = (TAp^rot- To es timate 5 TOt we have perfo rmed Monte 
Carlo simulations as in iNapolitano et al.l (|200ll ). be einnmg 
with a suite of analytical spherical stellar+DM models as 
described in N+0q 12 l, For each model with a fixed gravita- 
tional potential, we assume isotropy and an additional ro- 
tational component that increases with radius, then solve 
the Jeans equations and project to (Ja p . Finally we exam- 
ine the factor <S rot that relates the rotating and non-rotating 
"measurements" <7a p , finding this simple approximation: 

5 Iot « 1 + 0.05 — , (7) 

which is calculated for an aperture of 1 Rgg and turns out 
to be valid for a large range of galaxy massej^l We therefore 
apply this correction to the observed <ta p before matching to 

results almost unchanged. Following T+04, we have also checked 
that alternatively using fixed apertures of 1".6 and 2". 2, the me- 
dian Tdyn values are overestimated by 3^2% and lii%> we U 
within the typical uncertainty on each estimate and on the sam- 
ple's global scatter. 

11 We added an estimated black hole mass Mbh in the Jeans 
equations, using Mbh pred icted by the correlation oo — AfeH in 
iFerrarese fc Merrittl ll2000h . The result was to decrease Td yn by 
2±1%. 

12 The multicomponent model, in their § 3, includes 
a Hernquist (1990) stellar distributuion and an NFW 
l|Navarro. Frenk fc White] Il996l : iNavarro Frenk fc White! Il997ll 
spherical DM halo. 

13 This correction is much smaller than that found by PS94 be- 
cause these authors did not take into account the variations of 
rotation with radius, nor of measurements made within an aper- 
ture rather than along the major axis. 
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Figure 8. Ratio between rotation and dispersion of ETG sample, 
as a function of luminosity. Red diamonds, red triangles and blue 
circles are respectively slow rotator Es, fast rotator Es and SOs. 
The solid red and blue curves are functional fits to the binned 
data, for Es and SOs respectively, while the red dashed ones are 
for slow and fast rotator Es. Finally, the dashed grey line sets the 
boundary between slow and fast rotators. 

the models in step (vi) above. The correction increases the 
inferred Td yn values since rotation at ~ 7? c ff coupled with 
the (3 — assumption depresses the central ay for a given 
mass profile: a rotational component must be subtracted 
from the right hand side of Eq. ((3j). 

The trend with luminosity for the rotation correction 
in our galaxy sample is shown in Fig. [8] implying Tdyn cor- 
rections of ~ 1% for the brightest Es, and ~ 6% for the 
faintest SOs (as is well known, rotation is a stronger factor 
on average among fainter ETGs). When plotting results for 
all galaxies in the sample we use the median (Vniax/ '(to)-Lb 
trend to estimate their <5 rot . Where possible, we also classify 
the Es as fast- or slow-rotators, using Vmax/oo = 0.25 as 
the demarcation - a simple scheme that matches the more 
robust conclusions of C+06 in more than 90% of the over- 
lapping cases. 

4.2 Results from homologous luminosity profiles 

We can now derive Td yn , starting with the simplest case 
where we assume no rotation, and a homologous model for 
the luminosity distribution j*(r): an R 1 ' 4 profile which is 
completely determined by the known Lb and R e g for ev- 
ery galaxy. For the mass profile we assume initially the SIS 
model. Fitting the ao data, we show the mass-luminosity 
results in Fig. [5] (left panel). Binning the data, we fit the 
median relation Td yn oc is 7dyn and find 7d yn = 0.21 ±0.01. 
This is identical to the one for the E subsample alone, 
7d yla = 0.21 ± 0.01 (see Table [TJ. The SOs show a larger 
scatter and have a global slope of 7d yn = 0.18 ± 0.03, which 
steepens for faint galaxies (Mb > —20.5, 7d yn ~ 0.3) and 
appears to flatten or even decrease at higher luminosities 
(7dyn ^ 0). These slope results are scarcely changed by in- 
cluding the rotational correction (see Table [TJ, although the 
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Figure 9. Dynamical M/L in B-band within R c g as a function of luminosity assuming j* homology and an SIS total mass profile. Red 
squares and blue stars denote E and SO galaxies, respectively. Points with error bars are the median values and ±25% scatter for the 
galaxies in luminosity bins. Left panel: No rotation assumed. Linear best fits to the binned data are overplotted as straight lines. Right 
panel: Correction made for rotational support. 



Table 1. Slope of M/L-Lb relation for Es and SOs and different 
dynamical models (the Sersic profiles assume the n-Lg relation 
discussed in the text). The first five rows are the slope for the total 
(dynamical) mass, and the last row due to stars only, as derived 
with our stellar-populations model (Section [3). *The faintest SOs 
have a steeper slope than the brighter ones (see Fig. Uncer- 
tainties on slopes are the la scatter computed by a bootstrap 
method. 



Model 


" E 


7S0 


7tot 


R 1 / 4 +SIS 


0.21 ± 0.01 


0.18 ± 0.03 


0.21 ± 0.01 


ijl/4 + SIS+rot 


0.20 ± 0.01 


0.18 ± 0.03 


0.20 ± 0.01 


JJ^ + const-M /L + rot 


0.20 ± 0.01 


0.18 ± 0.03 


0.20 ± 0.01 


Scrsic+SIS+rot 


0.21 ± 0.01 


0.20 ± 0.03 


0.21 ± 0.01 


Sersic+const-.M/ L+rot 


0.10 ± 0.01 


0.10 ± 0.03 


0.13 ± 0.02 


Stars 


0.02 ± 0.01 


0.17 ± 0.03 * 


0.06 ± 0.01 



normalization of Td yn is increased for the SOs (see Fig. [9] 
right), an issue to which we will return in Sj5] 

Previous dynamical studies of ETGs using j* homology 
have found a v ariety of tilt slopes, ranging from 7dvn ~ 0-1 
to ~ 0.3 (e.g. IJorgensen et all Il993l ; iBernardi etafl 120031 ; 
IPadmanabhan et al.ll2004l ; T+04). The average of the lit- 
erature B-band values in Table 1 of PS96 yields 7d yn = 
0.25 ± 0.05, consistent with our result. 

Now considering the other extreme assumption for the 
mass model, const ant- M/L, the extrapolation to R c s after 
fitting to (To changes the Td yn normalization, correspond- 
ing to K = 2.05 at i? cff in Eq. (0 for SIS, and K = 1.93 
for constant- M/L. The slope of the Td yn -is relation is on 
the other hand unchanged (see Table [1]). The K difference 
does raise the interesting possibility of mass profile non- 
homology, e.g. a systematic change with luminosity. As with 
the IMF toy model in Section 13.31 we can consider an ar- 
bitrary case where the faintest galaxies have constant-M/L 
profiles, and the brightest ones have SIS. This would increase 
Yd y n by ~ 0.02, i.e. mass non-homology does not appear to 



be a significant contributor to the FP tilt, assuming j* ho- 
mology. 



4.3 Results from generalized luminosity profiles 

We next relax the j* homology assumption, allowing for 
more realistic luminosity profiles based on the Sersic law, 
with surface brightness profiles expressed as: 

Li(R)(xC-{R/R cH ) 1/n , (8) 

where C is a constant and n is an index of pro- 
file curvature which correlates with luminosity, such 
that the brighter gal axies have higher n (less curved 
profil e s; e.g. Caon. Capaccioli fc D'Onofriol 19931; Graham! 



I l998l; iGraham fc Guzman! 120031 ; iMamon fc Lokasi I2005I 
iKormendv et al.ll2008l ). As illustrated by Eq. ([5}, for a given 
dispersion profile, changing the shape of j, (r) will affect the 
inferred mass. Thus it is important to explore the impact 
of j, non-homology on Td yn , which may be expressed as a 
trend with luminosity K = K(n\Ls)- This is all a fancy 
way to say that accurate dynamical results require accurate 
luminosity profile models. 

The u-Lb correlation has been investigated for our 
galaxy sample by PS97. From the overall ETG sample in 
their fig. 5, we define a simple relation where n ~ L B 2 for 
M B > —20, and n — 4 for all the brighter galaxieiFr This 



14 PS97 noted that at least one other study found higher val- 
ues of n for the brightest galaxies, but commented that those 
results were more sensitive to the outer profiles than to the cen- 
tral regions of relevance here. Si milar concerns might ap ply to the 
recent smaller galaxy sample of Kor mendv et all (120081 ) , but it is 
beyond the scope of our paper to re-investigate n dependencies 
in detail. If n were systematically higher for the brighter galaxies, 
then these systems' T^yn results would be lower (cf. next foot- 
note) .Note also that the R e g values that we use were obtained by 
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relation also applies for the E subsample, and we assume 
that it does for the SOs as well. 

We now use the ti-Lb relation to construct the j*(r) 
Sersic profile for each galaxy as needed for the dynamical 
modeling (Section 14. Since we have examined the effects 
of rotation in Section T4. 21 we will here skip over the simpli- 
fied case of no rotation. The resulting Td yn values for both 
SIS and constant-M/L cases are summarized in Tableland 
Fig. If 01 For the SIS case, relaxing the homology slightly 
changes the slope 7d yn - However, in the constant-M/L case, 
both the luminosity and mass profiles are affected, and sig- 
nificant differences arise. The masses are increased for the 
fainter galaxief^l. causing the Td yn slope to become shal- 
lower overall (7d yn = 0.13 ± 0.02), and even constant at 
lower luminosities (7 dyn = 0.05 ± 0.04 and 0.23 ± 0.02 for 
the faint and bright Es, respectively). The Es and SOs are 
again not noticeably different in their region of luminosity 
overlap. In Appendix|B]we investigate systematic uncertain- 
ties in these results, whose impact we will consider in the 
next section. 

Now reviewing the results of this and the previous Sec- 
tions, with the j, homology assumption, the steep Td yn 
slope relative to Y„ (7d yn = 0.20 vs 7* = 0.06) would 
imply that ~ 75% of the FP tilt is related to DM con- 
tent or some other factor. Including the (realistic) j* non- 
homologies changes the picture somewhat: if all galaxies 
have SIS mass profiles, the previous conclusion is unchanged. 



ij 1 / 4 fitting in PS96 rather than the self-consistent Sersic values, 
which could in principle affect the results for the fainter galaxies. 
15 This is because lower n for the fainter galaxies implies shal- 
lower central profiles of both luminosity and mass density, and 
therefore lower central velocity dispersions at a fixed mass, fi- 
nally requiring higher model masses to match the observations. 
This effect might be somewhat reversed by the central cusps of 
light generally found in faint Es to be superimposed on their over- 
all Sersic profiles, in ma ny cases on spatial scal es comparable to 
the fro measurement (e.g. lKormendv et al.ll2008f) . but it is beyond 
the scope of this paper to consider this possibility in detail. 



If they have steeper mass profiles, then the dynamical con- 
tribution to the tilt decreases, and for the fainter galaxies 
may even disappear. 

Thus our results suggest overall that DM contributes to 
the tilt for the brightest galaxies, while the contribution for 
the fainter galaxies is unclear but probably less. This con- 
clusion differs from that of T+04, who found using similar 
Sersic models and assuming constant-M/L, no need for a 
correlation between DM fraction and luminosity. 

Their galaxy sample is fainter and much smaller, so 
their results are actually consistent with ours in general. 
The exception is for the brightest galaxies, where the higher 
n-values of T+04 lead to less tilt than we find. In any case, 
it should be noted that reproducing the FP tilt without DM 
variation is not a unique solution, and as we have shown, 
DM could still be a primary driver of the tilt. 



5 DARK MATTER FRACTIONS 

Having analyzed the trends for stellar and total mass in 
our galaxy sample, we now examine the implications for 
DM content. We define the DM fraction within the three- 
dimensional radius r — 1 R c g by: 



/dm = 



M tl 



M* 



= 1- 



T 



d y n 



(9) 



where for physically meaningful results we should have 
T» ^ Tdyn and thus /dm 0. Strictly speaking, our de- 
rived T* should be deprojected before computing /dm, but 
we do not have the information necessary to do so. Given 
the negative colour gradients in ETGs, we expect the de- 
projected T* to be somewhat higher than in projection, and 
thus the true /dm to be somewhat lower. For this reason and 
especially because of the large IMF uncertainty, the absolute 
values for /dm are not definitive, but instead the relative 
variations are more robust and are the focus of our study. 

We now consider the /dm trends found for our galaxy 
sample, taking as a default the T, estimates from the gener- 
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Figure 11. Trends of DM fraction with luminosity, using SIS and constant-M/L mass models (left and right panels, respectively). 
Symbols are as in Figs. l9 UlOI 




Figure 12. Trends of DM fraction with velocity dispersion, using SIS and constant- M/L mass models (left and right panels, respectively). 
Symbols are as in Figs. I MTuTfTTl 



alized (t ga i, t, Z) BC03-based stellar populations model, and 
the Tdyn estimates from the dynamical models using gener- 
alized luminosity profiles (Section l4.3|l . As shown in Fig. 111! 
/dm increases with luminosity in the E galaxy subsample, 
but is constant or even decreasing for the SOs. The combined 
ETG sample has /dm increasing overall, but with the hint 
of a slope change at Mb ~ —20.5, from roughly constant at 
faint magnitudes to steeply increasing for brighter objects; 
the trends for the SOs and the correspondingly fainter Es 
are roughly consistent. These conclusions are valid for both 
bracketing mass profile cases (SIS and constant-M/L), al- 
though the slope change is less apparent for the SIS model. 
To quantify this breakdown, we have measured the slopes 
7dm of /dm-Lb relation for the two models, and found 
that 7dm ~ 0.5 for log is J>10.4L Q and 7dm ~ for 
log Lb < 10.4 Lq, clearly inconsistent within the errors. 



The trends with luminosity are mirrored by similar cor- 
relations with the velocity dispersion as we show in Fig. 
1121 Here the slope is steeper in general because of the com- 
bined effect of the Y*-<ro correlation shown in Fig.[7]and the 
stronger dependence between Td yn and (To- 

The DM fraction is typically /dm ~ 0.3 assuming a 
Salpeter IMF, with a broad range for individual galaxies 
from ~ to ~ 0.9 (rms scatter of ~ 0.15). About 15% of 
the galaxies have, within the errors, /dm < (typically 
those with low surface brightness /u e ff), an unphysical result 
which may indicate that the Salpeter IMF is inaccurate (cf. 
C+06); adopting a Chabrier IMF would imply more DM, 
with /dm ~ 0.6 typically, and only a tiny handful of galaxies 
with /dm < 0. Changing the IMF also flattens slightly the 
luminosity dependence of /dm, since this quantity is not 
directly proportional to T*. 
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Figure 13. DM fraction trends for different galaxy types, using SIS and constant- M/L mass models (left and right panels, respectively). 
Dashed lines show median binned values and shaded areas show ±25% scatter. Light grey shows slow rotator Es, medium grey shows 
fast rotator Es and dark grey shows SOs. 



We next look for any DM differences between the fast- 
rotator and slow-rotator Es, following the classification in 
Section 14.11 However, as shown in Fig. 1131 there is no dis- 
cernible difference; slow and rotators typically have /dm ~ 
0.35 and ~ 0.25, respectively, but this is consistent with a 
simple luminosity effect, since fast rotators are fainter on av- 
erage than slow rotators. This result appears contrary to the 
finding of C+06 (based on more detailed dynamical models 
and somewhat different stellar populations constraints for a 
much smaller galaxy sample) that there is an /dm disconti- 
nuity between slow and fast rotators. 

We also compare SOs in Fig. 1131 where it appears that 
their declining trend of /dm with luminosity is inconsistent 
with the Es in the same luminosity range. However, we cau- 
tion that our spherical dynamical models are most question- 
able for the SOs, so the overall situation appears consistent 
with a continuous trend of /dm with luminosity for all ETGs, 
independent of morphology and rotation. For the rest of the 
paper, we will therefore generally lump all these ETG sub- 
classes together as one population. 

Before continuing further, we check once more the ef- 
fects of systematic uncertainties, as detailed in Appendix [Cl 
Despite the uncertainties, our default model is consistent 
with results on DM content at larger radii, and we therefore 
consider the overall mild increase of /dm to be robust, with 
the inflection at intermediate luminosities perhaps less so. 

How do our results compare to previous studies of DM 
trends in ETG centres? The analysis most similar to ours 
is from T+04. As discussed in Section 14.31 they found no 
indication of a correlation between /dm and luminosity, but 
their sample was primarily of faint galaxies, where we also 
found the correlation is weak. If we adopted higher Sersic 
indices for the brightest galaxies, the correlation would also 
weaken for them, but this scenario would seemingly be in- 
consistent with large- radius tracers of DM (see App. [Cjl . 

iBorriello et~aT1 (120031 ) modelled a large sample of ETGs 
dynamically and claimed that the flatness of the FP would 
not permit centrally-concentrated DM halos as predicted by 



cosmological models. Their results imply 7* = 0.27 ± 0.04, 
thereby explaining all the tilt through the stellar populations 
- in flat contradiction to our 7, = 0.06±0.01. This is mainly 
the consequence of their choices to not allow for a systematic 
variation of the virial DM fraction with luminosity in their 
model, and to use homologous j*(r) profiles which we have 

seen produc e misleading results. 

Finally, IPadmanabhan et alj (|2004j ) analyzed a large 
sample of SDSS ETGs, using a combination of stellar pop- 
ulations and dynamical models. Although the different red- 
shift ranges make comparisons not straightforward, their re- 
sults do appear roughly equivalent to ours, with 7, ~ 
and 7dyn ~ 0.17, and even a hint of a flattening of the /dm 
slope at lower luminosities. Note however that they used an 
inaccurate homologous j»(r) profile. 



6 IMPLICATIONS: DARK MATTER AND 
GALAXY FORMATION 

The trends we have seen for /dm as a function of lumi- 
nosity could provide fresh clues to galaxy formation. The 
most basic interpretation of central DM variations is that 
they reflect variations in total DM within the virial radius. 
Assuming that the Universal baryon fraction is roughly con- 
served from galaxy to galaxy, the implication is then that 
higher /dm means lower efficiencies of star formation esf- In 
this respect, the trends we find are qualitatively expected. 
Both observations and theory point to a universal U-shaped 
trend of esf (or equivalently virial M/L) wi th luminosity, 
and a peak efficiency at M* ~ 10 11 (e.g . Benso n et al 



|2000l ; lMarinoni fc Hudsonl2002l ; N+05; ivan den Bosch eta! 
120071) . 

Physically, the lowest-mass galaxies are least able to re- 
tain their primordial gas content long enough to form many 
stars, since their gravitational potential wells are not deep 
enough to prevent ejection from supernovae feedback. More 
massive galaxies are increasingly able to inhibit feedback 
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and form more stars, but at a certain mass scale, addi- 
tional processes kick in such as AGN feed back, inhibiting 
gas cooling and decreasing ssf again (e.g. ICattaneo et al, 
20061 ; IShankar et aill2006l ; iKavirai et al.ll2007l ; llbrtora et al. 
2009). Thus, the lowest-mass and the highest-mass galax- 
ies are the most DM-dominated. Our current galaxy sample 
does not extend faint enough to discern any U-shape, but the 
change we see in the /dm trend below scales of Mb ~ —20.5 
or M* ~ 10 11 Mq does coincide with the generically expected 
minimum of DM content. The consistencies of the trends for 
the ETG sub-types (SOs, fast-rotator Es, slow-rotator Es) 
suggest that the dominant driver of star formation is mass, 
not angular momentum. 

It is of course a stretch to draw firm conclusions about 
virial quantities based on data from scales <C R e ff- The cen- 
tral DM content that we are actually probing may be decou- 
pled from the overall DM content in several ways: the central 
DM density reflects the ambient density at the time of ini- 
tial halo collapse; the baryons could have interacted with 
the DM and changed its distribution; and the /dm quantity 
that we measure is somewhat dependent on the particular 
values of R c g for the stars rather than simply probing the 
DM properties. To allow for such effects, and to provide 
quantitative marks for comparison to cosmological theory, 
we now consider the properties of the DM alone, in terms of 
its average density within some small radius, (pdm}. 

In order to make critical comparisons with a literature 
study discussed below, we estimate (pdm) within 2 J? e ff, ex- 
trapolating our usual models outwards in radius. We present 
this result vs stellar mass for the whole ETG sample in 
Fig. 1141 using a Salpeter IMF and alternatively the SIS 
or const ant- M/L mass profile. Although these bracketing 
mass profiles gave similar results for /dm within 1 R c g, at 
2 R e s they start to diverge more, giving noticeably differ- 
ent results for (pdm). The less massive galaxies have in- 
creasingly dense DM halos, apparently reaching a plateau of 
(Pdm) ~ 0.05 M Q pc" 3 at masses below log M, / Mq ~ 1 1 . 

As a reality check, we comp are the results f or flat tened 
ETGs in the Coma cluster from iThomas et al. (2008, here- 
after T+08), who used detailed three-integral axisymmetric 
dynamical models of stellar kinematics to decompose the 
galaxies into their stellar and DM mass components. Their 
(Pdm) values from their NFW halo model match up remark- 
ably well with our SIS-based results. T+08 fitted their data 
with a logarithmic density-mass trend which would imply 
very high central (pdm) for the faintest galaxies. However, 
as we can see in the Figure, such conclusions would involve 
extrapolating outside the mass range covered by the data, 
and in fact the T+08 results do show some sign of the den- 
sity plateau at small masses which we find. 

Now we calculate predictions from ACDM cosmolog- 
ical models, a doptin g an NFW density profile, and the 
iBullock et all 1)200 ll ) mass-concentration relation, as dis- 
cussed in N+05. The final parameter in this model is the 
mass ratio between stars and DM within the virial radius, 
taking plausible values of alternatively M*/M v - lt = 0.1 or 
0.01 (corresponding to esf ~ 6% or ~ 60% , respectively, 
for a baryon fraction of 0.16; see N+05 and ISpergel et ail 
2003). The results for (pdm) are shown in Fig. 1141 where 
the model predictions are seen to be fairly consistent with 
the observations, including the bend in the (pdm) trend at 
similar galaxy masses. This bend is not caused by anything 
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Figure 14. Central DM density in ETGs, versus stellar mass. 
The density is averaged within a radius of 2 -R G ff , and is in units 
of M Q /kpc 3 ; a Salpeter IMF is assumed. The 25 th -75 th per- 
centiles are shown as darker and lighter blue for the Sersic+SIS 
and const— M/L models respectively. Dots are the T+08 galaxy 
sample with the best fit shown as purple solid line. Orange and 
green solid lines are the ACDM predictions using egp = 60% and 
6%, respectively. The black solid line is the model assuming a 
varying egp (see text for details). The dotted curves represent 
the range of model predictions due to the scatter in the i? e ff - 
M* relation. 



intrinsic to the DM itself, but by the radius adopted for 
measuring the density. As shown in Fig. [1] (right panel), 
the mass-_R c ff relation for ETGs has a bend, which proba- 
bly explains not only the density trends seen in this section 
but also the /dm results of the previous section. For bright 
galaxies, i? e fi increases rapidly with mass, probing quickly 
into regions contain more DM, at lower averaged densities. 
Fainter galaxies have less quickly varying _R c ff which thus 
tracks the slowly-varying DM scale radius more closely, so 
that the observable DM properties are roughly constant. 

In more detail, the esf = 60% theoretical case coin- 
cides roughly with our observational findings assuming a 
constant- M/L profile, and the esf = 6% case coincides 
nicely with our SIS-based findings. However, neither of these 
cases is plausible observationally or theoretically for the full 
range of galaxy masses. Our final case invokes a transi- 
tion from high esf for the faint galaxies to low esf for the 
bright galaxies, which is generically expected from various 
lines of evidence. More specifically, motivated by the find- 
ings of N+05 based on radially-extended dynam ical studies 
of ETGs (see also e.g.. iNapolitano et al.ll2009al ), we assume 
that the bright galaxies have esF decreasing steadily from 
60% to 16% in the mass range of logM, = 11-12 Mq; 
the faint galaxies have a constant esf = 90%. As shown in 
Fig. 1141 this model would be consistent with our SIS findings 
at the bright end, and with constant-M/L at the faint end. 
However, the NFW-based models themselves would have 
roughly SIS profiles for the entire range of luminosity, which 
means this set of model assumptions is not self-consistent. 
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It is beyond the scope of this paper to explore the possible 
combinations of DM parameters that would be fully consis- 
tent with the data, but we speculate that the low-luminosity 
objects have low-concentration DM haloes. Note that chang- 
ing the IMF to Chabrier would not significantly change these 
conclusions, since the data curves would shift up and to the 
left in Fig. M 



7 CONCLUSIONS 

The relative amounts of dark and luminous mass in ETGs is 
crucial information for understanding the internal structure 
of these systems and their formation mechanisms. In this 
paper we have analyzed both the stellar and dynamical M/L 
in the central regions of one of the largest homogeneous 
samples of local early-type galaxies, provided by PS96. 

We estimate the stellar content by accurate stellar pop- 
ulation synthesis models of several observed colours using 
the BC03 prescription. We measure dynamical masses us- 
ing the observed central velocity dispersion <ro and several 
simplyfing assumptions in the Jeans equations. 

We find that the stellar M/L, T*, has a shallow trend 
with luminosity with a slope ~ 0.06 for the whole ETG 
sample (with SOs showing a steeper trend than the Es: see 
Table [TJ. Dynamical M/L, on the other hand, have a slope 
for the M/L tx Lb 1 relation of is 0.21 ±0.01 when consider- 
ing ETGs as a (photometrical and kinematical) homologous 
galaxy family, i.e. fully consistent with results derived in 
local galaxies' B-band FP. 

For the non-homology case (i.e. assuming the Sersic pro- 
file for the light distribution and differential rotation within 
Red), we find that using the SIS model as the total mass dis- 
tribution does not much affect the M/L slope and thus not 
the FP tilt either. On the contrary, non-homology can ac- 
count for as much as ~ 40% if considering the constant-M/L 
model, and even more (up to 80%) for the faint systems. 

A further 30% (i.e. 0.06/0.21) is provided by the T, 
slope. The residual contribution to the M/L oc Lb 1 slope (~ 
70% for the SIS model and 30% for constant-Af/L) is mainly 
due to a variation with luminosity of their DM fraction. 

It must be stressed that this average budget of 7 con- 
tributions masks a more complicated distribution with lu- 
minosity. For instance, for the bright/massive galaxies (i.e. 
log Lb J>10.4 1/0 and log M, ^11.3 Mq) which have a quasi- 
R 1 / 4 profile and little or no rotation, the effect of the non- 
homologies is minimal and the slope of the M/L oc Lb 1 re- 
mains steeper than the faint systems where non-homologies 
can account for almost all the slope 7. This, obviously, re- 
lates to the trend of the DM fractions discussed in Section [5] 

Here we have seen that /dm is strongly varying with 
luminosity and mass. In particular, we observe a dichotomy 
in DM content of bright and faint Es: galaxy brighter than 
Mb ~ —20.5 and more massive than log M, ~ 11 — 11.3 Mq 
have an increasingly larger /dm while galaxies lying below 
these luminous and mass scales invert the trend, such that 
/dm is constant or marginally decreasing with luminosity 
and mass. When separating the E sample into "slow" and 
"fast" rotators it is evident that this two-fold trend is mainly 
found in the fast rotator systems (see Fig. 1 1 3p . These two 
kinematical varieties do not show large differences in their 
/dm properties. In particular, we do not find significant ev- 



idence for systematically lower /dm(J!ch) for the fast rota- 
tor variety (C+06), although with a large scatter one might 
make such a conclusion using a small statistical sample. The 
inclusion of the ellipticity and orbital anisotropy would in- 
crease the steepness of the faint/less massive sample, but 
would leave unaffected the bright/massive galaxy range, still 
maintaining the dichotomy (Fig. IC1I bottom right). As an 
alternative to a variable DM content, we have briefly ana- 
lyzed the effect of a change of IMF as a function of luminos- 
ity (see Fig. [B), which could also explain the FP tilt. 

The /dm dichotomy adds to other well known 
ETG correlations as found in the /i e — R e s rela- 
tion, FJ, size-luminosity (or size-mass) relations and in 
the correlations of Sersic index with both galaxy size 
and luminosity, as d i scussed in Section s [5] 13.21 and 01 
llCapaccioli et aT, 1992, Prugnicl & Simienlll997l . IShen et al.l 
120031 . iMatkovic fc Guzman! 120051 . etc.). Our results mir- 
ror the DM content in the outskirts of galaxies, where 
variations of virial M/L as a function of mass and lu- 
minosity have bee n found both in simulations and obser- 
vational analysis dBenson et al.| l2000. Ma rinoni fc Hudson! 
I2002L I van den Bosch et al.l 120071 ) . A similar dichotomy in 
DM content is not observed for SOs, which are gener- 
ally fainter and less massive than Es and are strongly af- 
fected by rotational support (influencing the normalization 
of Tdyn). They have a slightly higher DM fraction and 
show a monotonically decreasing trend with mass and lu- 
minosity, consistent with what i s known for spiral galaxies 
|Persic. Salucci fc Ashmanlll993l '). 

A continuity in DM content of galaxy as a function of 
amount of rotational support is possibly shown in Fig. 1131 
where we plot DM fractions as a function of luminosity for 
slow and fast rotators and lenticulars. 

Looking at the average central DM density, (pdm), 
we have found that this quantity has a fairly small scat- 
ter within the ETG sample. Albeit model dependent - the 
Sersic+SIS model providing (/9dm) which are 0.2 — 0.4 dex 
larger than the ones obtained with the const— M/L - the 
overall trend of the galaxy distribution decreases mono- 
tonically with the stellar mass and luminosity in good 
agreement with independent results obtained by Thomas 
et al. (2008) for ellipticals in the Coma cluster. Our 
larger statistical sample, though, has allowed us to dis- 
cern the presence of a "knee" in the distribution (around 
the usual mass/luminosity scale at log M* ~ 11 Mq and 
log Lb ~ 10.4 Lq) where the relation of the more mas- 
sive/luminous galaxies bends to a steeper slope than the 
one followed by the less massive/luminous systems. We have 
shown that this "knee" can be explained with the change of 
the slope in the R a a — M» relation at log M» ~ 11 Mq. 

As a robust estimator of the central DM density, (pdm) 
can be compared against the expected values for standard 
NFW profiles. The match found is broadly good, with the 
results obtained assuming the Sersic+SIS model favoring 
high dark-to-luminous mass ratios, i.e. lower star forma- 
tion efficiencies, while the constant-Af/L models fit lower 
M v i r /M* values, i.e. higher efficiencies. In order to match 
up with the picture where galaxies hav e star formation ef- 
ficiencies varying with t he stellar mass dBenson et alj|200d . 
iDekel fc Birnboirnl feooe'), we have shown that the DM den- 
sity characteristics should change with the mass with low 
mass systems being surrounded by more "cored" haloes 
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(well approximated by the constant-M/L models) and high 
mass systems by "cusped" haloes (here reproduced by the 
Sersic+SIS profile). 

This DM non-homology could be a possible expla- 
nation of the "anomalously" low halo concentration pa- 
rameters recently found modeling intermediate luminosity 
galaxies, compared to the giant ellipticals showing "regu- 
lar" concentration as exp ected from the ACDM simulations 
l|Napolitano et al ]l2009bh . In t his respect a model like the 
Einas to profile (lEinastdll965l. bu t see also [ Navarro et al.l 
|2004 2008, ICardone et al l 120051 . iGraham et alj 120061 ) or 
a phenomenological mo del including a wid e range of in- 
nermost density slopes l|Tortora et al.l f2007) provides suit- 
able working hypotheses to test on larger data sample with 
extended kinematics (e.g. Atlas3E ^1: or the PN.S Ellipti- 
cal Galaxy Survey: 



l2009al ; ICoccato et al 



Douglas et al 



JOUEli 

2009) 



2007:; iNapolitano et all 
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APPENDIX A: SYSTEMATIC EFFECTS IN 
THE STELLAR POPULATIONS MODELS 

Here we examine the role of systematic uncertainties in the 
stellar populations results, using different as sumptions and 
basis m odels (see, also, similar analy sis in iRettura et al 



(2006), iKannappan fc Gawiseri £2007) and IConrov et al 



(2008)). First we consider our default model based on BC03, 
using three different parameterizations for the SFH. In our 
reference model, r and Z (as well as t ga i) are free parame- 
ters fitted to each galaxy; a more simplified model has fixed 
t — 1 Gyr and Z = Zq corresponding to typical values for 
the whole sample; an even simpler SSP model has r = 
and Z — Zq. As shown in Fig. I All fleft panel), the distribu- 
tions of T* in all these models are fairly similar, except in 
the simplest case which shows a stronger tail to low values 
of T„. The impact of these differences is shown in Fig. IA2I 
(left panel), where it can be seen that overly restrictive mod- 
elling assumptions compensate with large variations in T* 
and thus steeper values for 7,. 

Next we compare basis model variations, starting with 
BC03 and M05. The different input stellar models and treat- 
ments of the thermally-pulsing asymptotic giant branch 
phase (TP-A GB) lead to different Y, predictions for the 
same colours l|Maraston et al.ll2006rF1 . Comparisons of some 
parameters derived from the two models with solar metallic- 
ity, given the same colour data, are made in Fig. IA2I (right 
panels). The inferred ages agree very well for older popu- 
lations, with M05 returning ages up to ~ 10% higher than 
BC03, while for younger populations, up to ~ 30% lower. 
The difference stems from M05 predicting V — R and V — I 
to be redder for young populations and bluer for old, while 
B — V is redder for all ages. The implications for T* are 



17 Recent preliminary updates of t he BC03 mode l s have included 
an im proved TP-AGB treatment iTBru zual 2007; Emini an et al.l 
2008). The colours and T* predictions are altered, particularly in 
the near-infrared, but not substantially for Z 0.4 Zq. These 
new models are more similar to M05 than BC03 but for old ages 
resemble BC03 jMcGrath et al.ll2007l l. 
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Figure Al. Distributions of recovered T^ s from PS96 galaxy sample, for different stellar populations models and assumptions. Top 
left: BC03+Salpeter model with different parameter assumptions: r = and Z = Zq (red); r = 0.75 Gyr and Z = Zq (blue); and t 
and Z free to vary (green). Top right: BC03 and M05 models compared, assuming r = and Z = Zq. The main panel shows a Salpeter 
IMF, with BC03 in blue and M05 in orange; the inset panel shows a lower-mass IMF, with BC03+Chabrier in red, and M05+Kroupa 
in green. Bottom left: BC03 and M05 models compared (solid and dashed lines, respectively), with Z free to vary. Bottom right: BC03 
and BJ01 models compared, assuming r = and Z = Zq. The solid black lines show BJ01, and blue shows BC03, in both cases with 
a Salpeter IMF. Lower-mass IMFs are also shown: red is BC03 with Chabrier IMF, short-dashed black is BJ01 with "scaled" Salpeter 
IMF, long-dashed black is BJ01 with "modified" Salpeter IMF, solid dark grey is BJ01 with Scalo IMF, dashed light grey is BJ01 
with "top-light" IMF slope, and dot-dashed light grey is BdJOl with "top-heavy" IMF slope (the latter two cases using the PEGASE 
prescription). 



that agreement is good for T* ;>6 Tq, while for lower val- 
ues the M05 predictions are smaller by up to a factor of two. 
The extension of the M05 results to smaller values of T» can 
also be seen in Fig. IA1I The final impact of these systematic 
differences on the trend of T* with luminosity is shown in 
Fig. IA2I (left panel): M05 yields a more steeply increasing 
trend. 

We finally examine the SSP models from BJOT BJ01 
predict a tight correlation between T, and galaxy colourful: 
using relations for the three colours B — V, B — R and V — I 
and minimizing a \ 2 function we determine the best fitted 
T* for the BdJOl pescriptions which are shown in Fig. IA1I 
(left panel). 

Assuming a Salpeter IMF, the use of different stellar 



8 This was obtained for spiral gala xies but it h as been proven 
to work for ETGs as well fBJOl. feell et alj|2003h 



prescriptions (BC03 vs BdJOl) has a negligible effect on 
the bulk of the T» distribution (e.g. solid black lines and 
blue ones in the Figure). Some of the assumed IM Fs in 
BdJO l ("scaled" and "modified" Salpeter and the IScalol 
1 19861 . for further details see BJ01) predict lower T*, 
with the Scaled Salpeter and Scalo IMFs giving similar 
results of the BC03+Chabrier one (red curves in the 
same Figure). Finally, u s ing B J01 results for PEGASE 
l|Fioc fc Rocca-Volmeragel Il997l ) prescription we obtain 
that: 1) a top-heavy IMF with a slope —0.85 gives T* 
values which are in the between of the Kroupa (2001) IMF 
(or Chabrier or Scalo) and Salpeter IMF predictions, while 
2) a top-light IMF with a slope —1.85 give much larger 
T, values. Note that distributions using directly BC03 
(red and blue lines in left panel of Fig. IA1|) have a larger 
spread around the peak distribution than the BdJOl results. 
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Figure A2. Left panel: Median binned stellar M/L as a function of luminosity, for our ETG sample and under various modelling 
assumptions. Different stellar populations prescriptions are shown (VZ96; BJ01; BC03; M05), as well as different parameter assumptions 
— see figure legend for details. Right panels: Comparisons of model parameters derived from the same data using different basis populations 
models. The top right panel shows the age, and the bottom right is stellar M/L, with the horizontal axes providing the results from BC03 
models, and the vertical axes M05 results — in both cases assuming a Salpeter IMF. The red circles are E galaxies, blue diamonds are 
SOs, and black stars are galaxies of other classifications not used in our main study. Solid lines show the one-to-one consistency relation. 
Grey squares in the bottom right panel show BJ0I modelling results. 



BJ01 results have been plotted in middle panel of Fig. 
I A2I like grey points. The slope of the relation shown in this 
figure is unchanged if we use the various prescriptions ana- 
lyzed in the paper above (see distributions using a Salpeter 
IMF in left panel of Fig. IA1|) ; on the contrary, a little change 
in the zero point is observed. 

As a final test, we compare results using different mod- 
els and data on the same galaxies. As a stellar synthesis 
model, C+06 fit single burs t models (using stellar prescrip- 
tion in IVazdekis et"ai1ll996l ). to some line-strength indices. 
Their estimates are on average 20% larger than ours (with 
a scat ter of 17%). Th i s disc repancy could not be fixed by 
fitting IVazdekis et~aH l|l996|) or BC03 SSP models to our 
galaxy colours. Some systematics can be ascribed, partially, 
to the extrapolation of line-strength indices (and velocity 
dispersion) from the very central regions to the effective 
radius, if some change in the average stellar population is 
present and unaccounted. 



APPENDIX B: INDEPENDENT CHECKS ON 
DYNAMICAL MASSES 

Given the simplifications of our Jeans models used to de- 
rive the dynamical masses (spherical quasi-isotropic mod- 
els), we test here using independent results whether our 
methods have introduced any systematic bias for T dyn - Our 
first cross-check is with C+06, who constructed detailed two- 
dimensional models of nearby ETGs. Our sample has 18 



galaxies m common with theirsH The main differences be- 
tween the two datasets are: 1) our distance moduli are on av- 
erage larger (0.05 mag) than C+06 but consistent within the 
scatter; 2) our effective radii are on average 5% larger with 
a median scatter of 16%; 3) the central velocity dispersions 
fro m C+06 are lower than the PS96 values by 6 ± 15 kms -1 
(see lEmsellem et al.ll200j ). 

C+06 constructed flattened, axisymmetric, constant- 
M/L dynamical model s, using both two-inte gral Jeans mod- 
els and three-integral ISchwarzschildl |l979) orbit models. 
Their luminosity models are multi-Gaussian expansions of 
the observed surface brightness profiles, and thus quite non- 
homologous. Converting our constant-M/L Sersic-based 
Tdyn results to the /-band and to the C+06 distances, we 
compare to their results in Fig. IB1I The masses are broadly 
consistent, with a systematic trend for ours to be higher by 
~ 20%. 

There are several possible reasons for this residual dis- 
crepancy, including rotation, orbital anisotropy variations, 
and galaxy flattening - all of which were handled in rigorous 
de tail by C+06 but not by our models. Based on the results 
of ICappellari et all l|2007h , the anisotropy effect should not 
correlate strongly with luminosity, but rotation and flatten- 
ing probably do. We also compare our modelled values of <r e 
(the velocity dispersion integrated over a 1 R e s aperture, 



19 Another 7 from their sample did not have measured B — L 
colours available for making a proper comparison between their 
/-band and our B-band Td yn results. 
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Figure Bl. Comparison of dynamical M/L estimates from our models (horizontal axes) and those from the literature (vertical axes). 
The diagonal lines show the one-to-one relation. Left panel: Comparison to C+06 results in the /-band. Open red and filled blue symbols 
are for three-integral Schwarzschild and two-integral Jeans dynamical models respectively. Right panel: Comparison to MD07 results in 
the B-band. Red symbols represent galaxies with results including C+06; blue circles show those based solely on also data from C+06, 
while those with blue circles use only results from C+06. 



folding in the rotational contribution) with their observed 
values, to see if our extrapolation from the central aperture 
could be generating the discrepancy. However, our a e values 
turn out to be lower by 13 jig kms -1 , which goes the wrong 
way to explain our higher masses. 



Next we consider the detailed spherical dynamical mod- 
els of G+01, with 16 galaxies in common. After shifting to 
the same distance scale, our Td yn values at 72 c ff are 27%±8% 
lower on average than theirs. Since their sample was focussed 
on round galaxies, we suspect again that flattening is play- 
ing a key role in the accuracy of our results, but that we 
have been able to largely compensate for its effects in our 
simplified modelling. 



finally we turn to th e dy namical results of 
Ivan der Marel fc van Dokkuml (|2007l . hereafter MD07) 
who com piled Td vn for 60 lo cal galaxies from t he lit- 
erature (Ivan der Marell 1 199 it iMagorrian et all 1 19981 : 



iKronawitter et alj l2000l : bebhardt et alj|2003t C+06). The 
original works made use of various types and quality of data 
and dynamical models, but should in general be superior to 
ours. The Td yn values are combined after homogenizing the 
distances and cosmology, and converting to the _B-band. As 
shown in Fig. IB 1 1 (right panel) , there is good agreement in 
general, but again a tendency for our results to overestimate 
the masses by ~ 20%, which appears to be an effect of 
the fainter, flatter galaxies. Note that our extrapolated 
a e shows no systematic offset from MD07. In order to 
potentially correct for a systematic error in our dynamical 
modelling, we adopt a heuristic correction of 66% and 90% 
for faint and bright galaxies separately. 



APPENDIX C: SYSTEMATIC 
UNCERTAINTIES FOR DARK MATTER 
FRACTION 

We consider finally how various systematic uncertainties 
could impact the /dm determinations. We first consider the 
stellar populations models. As detailed in Appendix [X] the 
model prescription that is used can have a noticeable effect 
on the T* trends. We show in Fig. IC1I (upper left panel) the 
differences engendered in /dm by adopting different models. 
Among the most plausible models, the results are roughly 
consistent, but the trend of /dm with luminosity tends to 
flatten or steepen with the use of M05 or BJ01, respectively, 
rather than BC03. 

We next consider uncertainties in the dynamical mod- 
els, starting with the assumed mass profile. As shown in 
Fig. IC1I (upper right panel), the bracketing models of SIS 
and const ant- M/L produce similar results for /dm- Testing 
the possibility that our n = 4 Sersic index for the bright Es 
is inaccur ate, we alternatively take the h i gher n — Mb rela- 
tion from lCaon. Capaccioli fc D'Onofrid (| 19931 ) as reported 
in PS97, and find that in a constant-M/L case, Td yn for 
the brighter galaxies decreases and the overall /dm trend 
is constant with luminosity (Fig. IC1I bottom left panel). 
However, an SIS profile is probably a better match for these 
galaxies, and in this case, changing n would not affect the 
results. Finally, we try to calibrate out the inaccuracies in 
our simplified Jeans modelling, based on the MD07 results, 
and find that the Tdyn values for the fainter galaxies might 
actually be lower, and the /dm slope with luminosity there- 
fore steeper (Fig. IC1I bottom right panel) . 

To quantify the effect of ellipticity (e) on our estimates, 
we have also selected E galaxies with e < 0.3 (as derived by 
RC3) . For these systems, the results are still consistent with 
an increasing (flat) trend of /dm with luminosity for bright 
(faint) galaxies. 

In summary, there are several potential competing sys- 
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Figure CI. Effects of systematic modelling uncertainties on DM fractions, for the overall EG sample. Unless otherwise stated, the mass 
model is constant-M/L, and the IMF is Salpeter. Top left panel: Changing the stellar populations basis model (see left panel in Fig. IA2| 
for line definitions). Top right panel: Changing the dynamical mass model from S IS (black) to constant-M/L (grey). B ottom left panel: 
Changing the Sersic index n in the dynamical modelling (black: PS97 values; erev: ICaon, Capaccioli & D'Onofrio (1993) values). Bottom 
right panel: Calibrating the dynamical models using MD07 (black: original; grey: recalibrated). Here a Chabrier IMF is used to avoid 
negative /dm values. 



tematic effects, and it is not clear which one might win out 
in biassing the /dm slope. Given this uncertainty, we carry 
out a different, critical test of confidence in our results. Find- 
ing results in the literature for the mass content of galaxies 
in our sample at large radii, we construct the M/L-gradient 
parameter VfT introduced by N+05. This simple but pow- 
erful metric is calculated from dynamical measurements of 
M/L at inner and outer radii by the following formula: 

V * T = /? Rc \ Hr* - 1) • ( C1 ) 

Given the longer lever arm, Vf T when available tells us with 
greater security whether or not an object is rich or poor in 
We compare /dm and V^T in Fig. IC2I and confirm 



that high- /dm objects from the current paper generally have 
high halo DM content in the literature while low- /dm have 
small T consistent with a lower global DM content. 



20 No attempt is made here to decompose the M/L measure- 
ments into stars and DM, i.e. to determine /dm- Instead, the 



broad premise is that T(jy n increases more rapidly with radius in 
galaxies with higher /dm- 
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Figure C2. Af/L-gradient parameter based on extended dy- 
namics, compared to central DM fraction (default SIS model 
with Chabrier IMF). Red and blue dots are E and SO galax- 
ies. Most of the data are taken from N+ 05, with several updates 
and additions from more recent literature IT rcodorcscu ct al. 2005; 
Schuberth et al.ll2006l; iDouglas et al.ll2007l; IWeiimans et alJhuOS ; 
Napolltamo'ljt"lrtTl20*09a! ~ lR^ 
20091) . 



